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ABSTRACT 


Recent data from SAS-2 on the galactic y-ray line flux 
as a function of longitude reveal a broad maximum in the 
y-ray intensity in the region |&|<30°. These data, when 
unfolded, imply that the low-energy (1-10 GeV) galactic 
cosmic-ray flux varies with galactocentric distance and is 
about an order of magnitude higher than the . local value in 
a toroidal region between 4 and 5 kpc from the galactic 
center. We further show that this enhancement can be 
plausibly accounted for by first-order Fermi acceleration, 
compression and trapping of cosmic-rays consistent with 
present ideas of galactic dynamics and galactic structure 
theory. Our calculations indicate that cosmic-rays in the 
4 to 5 kpc region are trapped an accelerated over a mean 
time of the order of a few million years or about 2 to 4 
times the assumed trapping time in the solar region of the 
galaxy. Cosmic-ray nucleons, cosmic-ray electrons and ion- 
ized hydrogen gas are found to have a strikingly similar 
distribution in the galaxy according to both the observa- 
tional data and the theoretical model discussed here. 
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INTRODUCTION : 

Recent observations of galactic y-radiation above 30 
MeV obtained from the second Small Astronomy Satellite (SAS-2) 
as reported by Kniffen, et al. (1973) have confirmed the 
earlier results of Kraushaar, et al. (1972) indicating a 
large increase of y-ray flux in the direction of the gal- 
actic center. The data obtained from SAS-2 with somewhat 
better angular resolution then that obtained by the OSO-3 
detector, revealed a broad , fairly flat maximum in the in- 
tensity distribution of the radiation from the galactic plane 
extending over a 60° longitude interval centered on the gal- 
actic center (i.e., j&j < 30° ) J ' . It has been shown previous- 
ly that the flux of radiation from the galactic plane in the 
disk region away from the galactic center ( | Jl | >, 60° ) can 
be accounted for as arising from the interactions of cosmic 
rays having the intensity observed in the solar region of the 
galaxy with interstellar gas having the column density of neu- 
tral hydrogen observed in the galactic plane (Kraushaar. et 
al „ 1972, Stecker 1973 (and references therein), Kniffen, 
et al. 1973). In our previous letter (Stecker, et al. 1974), 
we showed that the y-ray spectrum in the inner galaxy as 
reported by Kniffen, et al . (1973) could be described as 

^We will adopt here negative longitude designations such as 
-30° instead of the more conventional 330° because such de- 
signations lend themselves more naturally to a mathematical 
discussion of the flux distribution owing to a large amount 

of symmetry around the galactic center. 
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consisting of 70 per cent tT°-decay radiation produced from 
the same type of interaction process which could account for 
the general disk radiation with an additional 30 per cent 
contribution that can be accounted for as arising from 
Compton interactions of cosmic-ray electrons with the larger 
radiation fields existing in the inner galaxy. In that 
letter, an explanation was proposed for the enhanced y-ray 
flux from the inner galaxy as primarily due to an enhance- 
ment in the cosmic-ray intensity in the inner galaxy caused 
by the large-scale acceleration and compression of cosmic- 
rays. We will undertake here a more detailed examination of 
that picture based on an actual unfolding of the SAS-2 
longitude distribution data. 

II. UNFOLDING OF THE LONGITUDE DISTRIBUTION OF THE GALACTIC 

RADIATION AND THE IMPLIED COSMIC-RAY DISTRIBUTION. 

If we assume that the y-ray production rate in the 
galactic plane Q is only a function of galactocentric dis- 
tance W, we have basically a two-dimensional problem in 
cylindrical coordinates with no azimuthal dependence and 
symmetry in the longitude distribution about 0° longitude. 
The SAS-2 data are suggestive of some asymmetry (Kniffen, et 
al. 1973, Stecker, et al. 1974) but such asymmetry does not 
appear to be pronounced. With the assumption Q = Q (tff) , the 
y-ray line intensity in photons per cm * s ’ rad as a 
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function of galactic longitude can be expressed in the 

form 


i y W = 


2tt 


j. 

/ 


M 


db 


/ 


OVRq) cot b 
dp Q(ST) 


(1) 


where R © is the galactocentric distance of the sun (10 kpc) 
and h is the half-width of the gas disk of the galaxy. The 
quantities to (galactocentric distance) and p (heliocentric 
distance) will be given in natural galactic units (1 galactic 
unit = 10 kpc) so that p = cos it± (W^ - sin 2 £) . The 

galactic latitude is denoted by b and we assume that the 
observations are made with a detector which integrates the 
flux over all contributing latitudes up to b^. 


For the purposes of analysis, we will assume that the 
production rate Q in the outer part of the galaxy is 
roughly constant and equal to Q q for 0.8 < S' < 1.5, neglec- 
ting some possible effects from the outer arms, and is 
equal to 0 outside of this toroidal region. In the inner 
galactic region, we will take Q to be equal to some 
function Q. to be determined for GKS&L, = 0.8 and Q. = 0 for 
^>(3^. The observed flux, I^(£) can then be written as the 
sum of an inner-galactic component 1^ and an outer-galactic 
component I Q where 


= *i + *o 


U) 


= V 27r / db / 


(h/R Q ) cot b 


dp Q o (W) 


( 2 ) 
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and may be expressed in the form 


I.(A) 


to 2 

h cos £ /* M o Ch (w) 

/ duT y 1 - 9 - 5 7 - 

2ir J 0 (l-QT) W -sinl) l/ 2 

sin^i 


(3) 


By expressing the observed flux in the form given by 
equations ( 2 ) and (3) , one can obtain a unique solution for 
the production rate in the inner galaxy, Q^(5y) which can be 
found by inverting the convolution integral in equation (3) . 
Equation (3) is a type of Abel equation which may be solved 
by the method of Laplace transforms. We then obtain the re- 
sult 



where n 


sxn 


a. 


In order to determine Q^, Q q is first chosen to fit the 
data obtained by Kniffen, et al. near the anticenter direction. 
Then using the interpolation of the SAS-2 data shown in 
Figure 1, 1^ is obtained from equation (2) . The production 
rate in the inner galaxy Q^(a?)is then calculated from equa- 
tion (4) using average values for the quantity [d/dri (-I^/cosA) ] 
averaged over 10 ° longitude intervals consistent with the 
form of the SAS-2 results. 1^ and the average value of 



- 6 - 


dl./d.^ are shown in Figure 2 and the results obtained for 
hQ(GT) are shown in Figure 3. We then define the quantity 
taken to be equal to the y-ray production rate per HI 
atom per second from cosmic-ray-gas interactions in the 
solar galactic region. This quantity, q Q is found to be 
equal to (1.3±0.2) x 10~ 25 s -1 from an analysis of recent 
accelerator data (Stecker 1973) . The cosmic-ray enhancement 
factor in the galaxy I cr /I Q (taken relative to the value in 
the solar vicinity) is then given by 

I/I = (hQ) (5) 

cr 0 

g © n HI 

where n„ T is the density of atomic hydrogen seen in 21 cm 
radio emission and £ is a factor representing the amount of 
cool atomic hydrogen and molecular hydrogen not seen in the 
21 cm radio surveys. For further analysis, we will assume 
that the value of (h£) is independent of galactocentric dis- 
tance GT. The distribution n HI (Gr) has been discussed by 
various authors (Kerr (1969) , Westerhout (1970) , Shane (1972) , 
Sanders and Wrixon (1973) . These data are summarized in 
Figure 4. Using the data given in Figures 3 and 4, we have 
used equation (5) to obtain the cosmic-ray enhancement 

factor (I /I ) as a function of galactocentric distance as 
C IT Q 

shown in Figure 5. The normalization I /i Q (£=*)=1 gives a 
value for h£ of - 140 pc for the intensity of y-rays observed 
by SAS-2. This is a reasonable value for the galaxy as 
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previously discussed by Stecher and Stecker (1970) . 

The values for hQ given in Figure 3 are evaluated from 
equation (4) by integrating inward toward the galactic 
center. For this reason, the relative uncertainty in the 
value of hQ increases near the galactic center as shown in 
Figure 3. Keeping this in mind, we can draw the following 
conclusions : 

1. The production rate of galactic y-radiation has a 
strong variation with galactocentric distance. Since the 
spectrum of this radiation indicates a cosmic-ray production 
origin, this implies that the cosmic-ray intensity is not 
uniform in the galaxy. 

2. The 1-10 GeV cosmic-ray intensity in the toroidal re- 
gion defined by 0.4<W<0.5 exhibits a strong maximum as 
shown in Figure 5, indicating the cosmic-ray intensity in 
this region is roughly an order of magnitude higher than in 
the galactic region in the solar vicinity. 

III. COMPARISON OF THE COSMIC-RAY AND y-RAY DATA WITH OTHER 
LARGE-SCALE FEATURES IN THE GALAXY. 

The radial distribution of y-radiation (Figure 3) and 
implied distribution of cosmic-ray nuclei (Figure 5) bear a 
striking resemblence to the distribution of another galactic 
parameter which is of importance to galactic structure 
studies, viz., the ratio of ionized to neutral hydrogen in 



- 8 - 


the qalaxy n As we pointed out in our previous 

Hi! AIJ - 

letter (Stecker, et al. 1974) , the overall fraction of hy- 
drogen gas ionized in the galaxy in the so-called "giant 
HII regions" has a similar radial distribution (Reifenstein 
et al. 1969, Mezger 1970) which has been related to the 
formation of young hot stars in spiral arms triggered by 
density waves in the density wave theory of spiral structure 
(see e.g. Oort 1973, Shu 1973). However, an increase in the 
low-energy sub-relativistic cosmic-ray intensity related to 
the increase in the 1-10 GeV cosmic-ray intensity deduced 
here may also play a significant role in accounting for the 
HII/ HI distribution as suggested to us by G. Field {private 
communication) . The correlation between I /I^ and HII/HI 
is quite dramatic as shown in Figure 6. The dashed line in 
Figure 6 also shows a function given by Shu (1973) to corre- 
late HII/HI with a compression factor given by the density 
wave model. The asymptotic dashed line between 3 and 4 kpc 
corresponds to the inner Lindblad resonance in this model. 

We have previously suggested that compression and accelera- 
tion can account for the enhancement of cosmic-rays implied 
for the 4 to 5 kpc region (Stecker, et al. 1974) as related 
to the "expanding arm" model (Oort 1970, Van der Kruit 1971) 
and will further elaborate on that hypothesis here, but we 
note that whatever the actual cause, be it expanding arm or 
density wave (Roberts 1970, 1973) it is the compression it- 
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self which can account for the increase in I /I~. 

There is a similar increase in observed non-thermal 
radio emission from high energy electrons toward the galactic 
center, but the situation with the non-thermal emission is 
more complex. Here also compression and acceleration effects 
acting on the cosmic-ray electrons together with an enhanced 
magnetic field strength are expected for the compressed 
region between 4 and 5 kpc and can be invoked to explain 
the data on non-thermal radio emission. 

In a recent review, Price (1973) has discussed the 
longitude distribution of the non-thermal radio background 
at 150 MHz. This distribution, as given by Price, is shown 
in Figure 7 along with the SAS-2 y-ray data and a model 
constructed by Price to account for the distribution of the 
radio data. One can see from Figure 7 that there is a gen- 
eral correlation between the y-ray and radio distributions 
for 1 1 1^30° which is not surprizing if the cosmic-ray 

electron and proton components are roughly proportional and 

. . 2 

the square of the magnetic field strength B ~n H * The corre- 
lation can be extended to the entire longitude range, if it 
is assumed that the peak in the radio intensity at i=0° can 
be attributed to a purely non-thermal radio emission region 
at the galactic center (Price 1973) . The rest of the 
emission has been accounted for by Price's model shown in 
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Figure 7 which supposes a general e flux distribution 

(similar to that we proposed for the y-radiation (Stecker 
et al. 1974)) with superimposed, relatively small variations 
due to smaller enhancements in the spiral arms. The maxima 
at 20 | <30 ° evident in the y-ray data are not obvious 

in the radio data. One reason for this may be the relative 
thickness of the radio disk of the galaxy as compared with 
the gas disk. 

THE ACCELERATION MODEL; 

An enhancement of cosmic-rays in a localized region 
of the galaxy can be due to three factors (1) an increase 
in the density of cosmic-ray sources (or alternatively the 
production rate) in the region (2) an increase in the 
trapping time (escape time) of cosmic-rays in the region, 
and (3) acceleration and compression of cosmic-rays in the 
region. We may write this as follows; 

£n(I cr /l t3 ) = + £n(T cr /T 0 ) + 6 (6) 

where the first term on the right hand side of equation (6) 
represents the enhancement in the production rate, the 
second term represents the trapping factor and the third 
time 6 represents the effect of acceleration and compression. 
The first term can only be guessed at, given our present 
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lack of knowledge of the ultimate origin of cosmic rays. 
However, as an indication of the possible effect of source 
density, one can note the distribution of supernova remnants 
in the galaxy which may be proportional to the density of 
cosmic-ray sources if we assume that cosmic-rays are 
produced by supernova explosions (Ginzburg and Syrovatsky 
1964) or in remnant pulsars (Gunn and Ostriker 1969) . 
According to Ilovaisky and Lequeux (1972) the density of 
supernova remnants increases roughly by a factor of 2 over 
the local galactic value for GJ^O.S and drops sharply for 
©>1.2. A more recent study by Clark, et al. (1973) is 
probably less susceptible to selection effects because 
their survey included remnants down to lower luminosity 
levels. The results obtained by Clark, et al . confirm the 
earlier results of Ilovaisky and Legueux regarding the 
general distribution of supernova remnants in the galaxy. 

We will therefore estimate here that Q. cr /Q© ~ 2 for the 
enhanced region 0.4<G><0.6. This factor, by itself, cannot 
account for the order-of -magnitude enhancement deduced for 

i. cr /%> ♦ 

The second factor, T /T^ is so difficult to estimate 
that we will treat it as a free parameter >1 to be solved 
for. It is not unreasonable to expect more effective 
trapping in the inner galaxy (i.e. T /T >1) due to com- 
pression and a resultant stronger magnetic field strength. 
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The factor 6 can be broken up into two parts, an acce- 
leration factor S ACC which accounts for acceleration of the 

more numerous lower energy cosmic-rays (given an assumed 

_r 

power-law differential energy spectrum of the form «E ) 
to an energy above the threshold for tt° production, and a 
simple density enhancement (<$ due to the lower volume of the 
compressed region in which the trapped particles find them- 
selves. If we designate the specific volume compression 
rate by 

a = -<1/V) (dV/dt) (7) 


then obviously 


S D - 


The acceleration factor can be estimated thermodynam- 
ically. Regardless of the exact details of this process, 
which may be a coherent, first order Fermi acceleration which 
can transfer the momentum of moving "clouds" with trapped 
magnetic irregularities to cosmic-rays trapped in the 
compressed region (Fermi 1954, Stecker et al. 1974), the 
increase in the energy of the individual cosmic-rays can 
be treated as an adiabatic compression heating of a "cosmic- 
ray gas molecule". The energy -volume relation then gives 
the energy enhancement factor as 


An(E/E 0 ) = (y-l)aT cr 


( 9 ) 
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4 5 

where y =j for relativistic cosmic-rays Y = y for sub- 

relativistic cosmic-rays. For a cosmic-ray differential 

_r 

energy spectrum of the form I “ E , we then find 


6 = (r-l)(Y-l)aT 

acc 1 cr 


and therefore 


(10) 


6 - fj(r-l) (y-l)+l]aT cr 

(Stecker, et al. 1974) . The value of the specific volume 
compression rate a(C3^ can be given in terms of the radial 
expansion velocity v r of the gas deduced from the 21 cm 
observations of the "expanding arm" feature. Under this 
interpretation 


R Q a(c3) = 


dv 


r - 


v 


dw 


r 

To 


The function of a(££T) deduced from the observations of v r (cr) 
tshane 1972, Sanders and Wrixon 1972, 1973) is shown in 
Figure 8. It is positive and maximal in the region of 
observed maximum y-xay emission (see Figure 3) and is nega- 
tive in the inner region uf< 3 kpc where there may be a 
significant drop in y-xa.y emission, although the uncertainties 
in that region are very large. 


Because of the large uncertainties in T and T ~ we 

cr vy 

will consider a range of values for T and T /T/-> . We 

cr cr u 

6 7 

will consider 3x10 <T cr <3xlO yr where the upper limit on 

T is taken to be the deduced age of the expanding feature 
cr 

(Oort 1970 ) , Van der Kruit 1971) . The cosmic-ray lifetime 
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in the local region, T 0 , is given by O’Dell, et al. (1973) 

6 7 

and Brown et al. (1973) , to be 10 yr<T <10 yr. Taking 

a=2xl0 - '*' 5 s -1 in the region of maximum compression (see 

Figure 8) and assuming I cr / I 0 s 15 and Q cr /£ 0 ~ 2 ' F-2.5, 

6“1.75<*T =1 . 5xlO - ^T , the value given for the parameters 

C c 

T /T^ , e and T and T consistent with equation (6) are 
cr' o' cr o 

6 

given in Table 1. The factor e represents the enhancement 
effect due to acceleration and compression. 


CONCLUSION 5 

Our results indicate that the galactic y-ray production 
rate, and by implication, the cosmic-ray nucleon intensity 
in the galaxy rises to a sharp maximum in a toroidal region 
between 4 and 5 kpc from the galactic center. This distri- 
bution is strikingly similar to that observed for ionized 
hydrogen in the galaxy and the radio-emission model of 
Price (1973) . However, no such correlation is observed 
with the neutral hydrogen distribution. The fact that the 
average neutral hydrogen distribution in the galaxy peaks 
in the outer region between 13 and 14 kpc from the galactic 
center and falls slowly inward toward the galactic center 
(see Figure 4) strengthens the argument that the cosmic-rays 
themselves must increase sharply in the 4 to 5 kpc region 
in order to explain the observed y-ray longitude distribu- 
tion. One should note that the strong 4 to 5 kpc maximum 
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appears to lie on the inside of an arm as is the case with 
the HI I maximum. 

It seems reasonable to expect relatively small pertur- 
bations superimposed on this general pattern due to some 
concentration of cosmic-rays in the spiral arms as in the 
model of Price (1973) . These perturbations have been neglected 
here along with possible fluctuations due to dense clouds 
(Black and Fazio 1973) and perturbations at galactocentric 
distances greater than 8 kpc. Therefore, our assumptions 
of constant h£ throughout and constant hQ for G3>0.8 should 
not be taken too seriously but they should represent average 
values indicating no strong features in this range. More 
detailed studies of galactic y-rays with better angular 
resolution should reveal the finer scale perturbations 
mentioned above. 

It thus appears that, as has been hoped by many over the 
last decade , y-ray observations have begun to reveal impor- 
tant information on the large-scale distribution of galactic 
cosmic-rays. This information would appear to support the 
galactic origin hypothesis for most cosmic-rays (see e.g. 
Ginzburg and Syrovatsky 1964) . 

Moreover, our calculations indicate that large-scale 
compression and first order Fermi acceleration may be 
occurring in the inner galaxy as previously suggested 
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(Stecker , et al . 1974). Such a model supplements our 
ideas on galactic structure as revealed by radio and optical 
data . 
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Table I. Lifetime of Cosmic-Rays in the 4 to 5 kpc Region 

Calculated from the Acceleration-Compression Model 


T cr/ T Q 

6 

e 

T cr (1 ° 6y) 

T (10 6 y) 

2 

3.7 

12.6 

6.3 

3 

in 

« 

CM 

8.7 

2.9 

4 

1.9 

6.0 

1.5 
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FIGURE CAPTIONS : 

Figure 1. Distribution in galactic longitude of the galactic 
y-ray line flux as observed by 0S0-3 (Kraushaar, 
et al. 1972) and SAS-2 (Kniffen, et al . 1973) 
normalized arbitrarily for purposes of comparison. 
Statistical uncertainties quoted for these results 
(not plotted) are significantly greater for the 
0S0-3 results as compared with the SAS-2 results. 

The interpolation of the SAS-2 data was obtained 
on the assumption of a smooth average variation 
in the y-ray production rate as discussed in the 
text. 

Figure 2. The function I^(&) as defined in the text obtained 
by using the SAS-2 data for £<50° and defining 
I.(A)=0 for £>60°. The values for Al./dl used 
for calculating the y- ray production rate 

(ft) were obtained from the slopes of the jagged 
line shown in the figure. 

Figure 3. The value for hQ(£0 calculated from equations (2) 
and (4) as discussed in the text. A constant 
value for hQ is assumed for Ro^Skpc. Uncertainties 
in the determination of the value of hQ ( ccr) as 
evaluated from the integral equation (4) grow quite 
large in the central region W<4 kpc . 
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Figure 4 . The neutral hydrogen density n (O') based on the 
discussion of Kerr (1969 ) , Shane (1972) , and 
Sanders and Wrixon (1973) . 


Figure 5. The calculated cosmic-ray distribution relative 
to the flux in the solar vicinity (I cr /I G ) . A 
mean constant value for R SJ> 8 kpc is assumed. 

For purposes of calculation, it is assumed that 
h£ is a constant which is evaluated to be -140 pc. 

Figure 6. The correlation between I cr (w) and HII/HI (w) is 
shown here. Also shown is a theoretical function 
used by Shu (1973) in discussing the observations 
of HII/HI as correlated with compression effects 
in the density wave theory of galactic structure. 

Figure 7. A comparison of the longitude distribution of 
galactic y-radiation and 150 MH& synchrotron 
radiation. The central peak in the synchrotron 
radiation distribution appears to originate in 
the central region of the galaxy. The remaining 
radiation has been interpreted in a phenomenological 
model given by Price (1973) (shown by the dashed 
line) which also appears to provide a good 
representation of the y-ray distribution. 

The specific compression rate a (IS) obtained from 
the radial velocity data v r (uJ) based on the obser- 


Figure 8. 
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vations of Shane (1972) and Sanders and Wrixon 
(1972, 1973) . 
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